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Abstract. This article addresses a problem of micromagnetics: the reversal 
of magnetic moments in layered spring magnets. A one-dimensional model 
is used of a film consisting of several atomic layers of a soft material on top 
of several atomic layers of a hard material. Each atomic layer is taken to be 
uniformly magnetized, and spatial inhomogeneities within an atomic layer are 
neglected. The state of such a system is described by a chain of magnetic 
spin vectors. Each spin vector behaves like a spinning top driven locally by 
the effective magnetic field and subject to damping (Landau-Lifshitz— Gilbert 
equation). A numerical integration scheme for the LLG equation is presented 
that is unconditionally stable and preserves the magnitude of the magnetiza- 
tion vector at all times. The results of numerical investigations for a bilayer 
in a rotating in-plane magnetic field show hysteresis with a basic period of 2ir 
at moderate fields and hysteresis with a basic period of 7r at strong fields. 

1. Introduction. Exchange-spring coupled magnets (spring magnets, for short) 
hold significant promise for applications in information recording and storage de- 
vices. Spring magnets consist of nanodispersed hard and soft magnetic phases that 
are coupled at the interfaces. (In a hard material, the magnetic moment tends to 
be aligned with the easy axis; in a soft material, it is more or less free to align 
itself with the local magnetic field.) The superior magnetic properties of a spring 
magnet stem from the fact that the soft phase enhances the magnetization of the 
composite jy, |[ g, [|, ||, ||] . Since the performance of a spring magnet is determined 
by the stability of the soft phase against magnetization reversal, it is important to 
identify the factors affecting the reversal process. 

Thin films provide an interesting class of simple models for which one can perform 
both physical and computational experiments. A spring- magnet structure can be 
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realized by interleaving hard and soft magnetic layers, and because the magnetic 
properties vary predominantly along the normal direction, the structure of such 
spring magnets is essentially one dimensional. 

In this article we investigate magnetic reversal in a hard/soft bilayer — a layer 
of soft material on top of a layer of hard material — with strong coupling at the 
interface. The hard and soft layers both consist of several atomic layers; each 
atomic layer is treated as uniformly magnetized, and spatial inhomogeneities within 
an atomic layer are neglected. The state of the bilayer is thus described by a chain 
of spins, each spin representing the magnetic moment of an atomic layer. 

The dynamics of a magnetic moment are entirely local. A magnetic moment is 
like a spinning top, which is driven by the effective magnetic field and subject to 
damping. The relevant equation was first formulated by Landau and Lifshitz 
and later given in an equivalent form by Gilbert ||. The local effective field is 
derived variationally from an energy functional ||. 

The Landau-Lifshitz-Gilbert equation preserves the magnitude of the magnetic 
moment, and it is important to maintain this property in numerical approximations. 
Our first purpose in this article is to present an integration scheme that conserves 
magnetization at all times. 

Our second purpose is to present some results of numerical simulations for a 
bilayer. They show two types of rotational hysteresis: one at moderate fields with a 
basic period of 2ir, which is associated with the irreversible behavior of the chirality 
("handedness") of the chain of magnetic spins in the soft layers; another at strong 
fields with a basic period of 7r, which is associated with the irreversible behavior 
of a full-length transition of the chain of magnetic spins in both the hard and 
the soft layers. These results explain the experimental observation of hysteresis in 
some torque measurements [ [To| . They also agree qualitatively with some magneto- 
optical measurements of the magnetization angle However, they differ at the 
quantitative level. The discrepancy is due to the one-dimensional model, which 
does not allow for the nucleation and motion of nanodomains. 

The remainder of this article consists of three sections. In Section ^, we describe 
the mathematical model, together with the numerical approximation procedure. In 
Section^, we present some simulation results for hard/soft bilayers. (Additional 
results are presented in fll2||.) In Section^, we summarize our conclusions. 

2. Mathematical model. A layered spring magnet is a multilayer structure, 
which consists of Nh atomic layers of a hard magnetic material adjacent to N s 
atomic layers of a soft magnetic material, 

Hard layers : i € Ih = {1, • • • , Nh}, 
Soft layers : i £ I 8 = {N h + 1, . . . ,N h + N s }. 

We put J = 4 U I s and N — + N s . The atomic layers are homogeneous, 
and variations occur only in the direction normal to the layers. We assume for 
convenience that the atomic layers are equally thick; their thickness d is of the 
order of angstroms (1 A = 1.0 • 1CP 8 cm). 

We adopt a right-handed Cartesian (x, y, z) coordinate system, where the x and 
y axes are in the plane of an atomic layer, the x axis coincides with the easy axis 
of the hard material, and the z axis is in the direction normal to the layers; e x , e y , 
and e z are the unit vectors in the direction of increasing x, y, and z, respectively. 
In a polar (</>, ff) coordinate system, (j> is the out-of-plane angle and 9 the in-plane 
angle measured counterclockwise from the positive x axis. 
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The state of the bilayer is completely described by the set of magnetic moments, 

M = {Mi-.ie /}. (l) 

Each Mi is a vector- valued function of time t, with components M,- lX , Mi tV , and 
Mi >z . The magnitude Mi of Mi is the magnetization, the unit vector m; = Mi/ 'Mi 
is the magnetic spin in the zth layer. The magnetization is constant at all times 
and equal to the local saturation magnetization, 



Mi(t) = Miiriiit), with Mi 



M s 



if i e 4, 
if i e I s . 



(2) 



Here, Mh and M s are the values of the saturation magnetization for the hard and 
soft material, respectively. Each magnetic spin can be specified in terms of its 
Cartesian or polar components, 



mi 



{m ltX ,m i ^,m ltZ ) t = (cos<^ cos( 



, cos 0i sin t 



, sine 



0* 



(3) 



Thus, 9i is the in-plane angle of m; with the easy axis of the hard material (mea- 
sured from the positive x direction), <j)i the out- of -plane angle of m,. 

2.1. Dynamics of the magnetic moment. A magnetic moment is like a spinning 
top, which is driven by the local effective magnetic field and subject to damping. 
The equation of motion for Mi is the Landau-Lifshitz-Gilbert (LLG) equation, 

3M, 



dt 



i e I. 



(4) 



Here, Hi is the effective magnetic field in the ith layer, 7 the gyromagnetic constant, 
and g a (dimensionless) damping coefficient. Note that the LLG equation yields a 
magnetic moment whose magnitude is constant in time. An equivalent form of the 
LLG equation is 

' >lVl ' -c (Mi x Hi) + -i-Mi x (Mi x H { ) 



dt 



Mi 



i e 1, 



(•5) 



where c = 7/(1 + g 2 ). We rescale t by a factor c and take c = 1 from here on. 

The effective magnetic field is found by taking the variational derivative of the 
free energy, 

H * = (6) 

where H a is the externally applied field, which we take to be uniform and constant 
in time. The free-energy density F is the sum of the exchange energy, the anisotropy 
energy, and the demagnetization energy, 



T\M] 



\Mz) 



dm 



dz 



K(z) \m x e x \ 2 + ~(4tt)(M ■ e z 



(7) 



Here, Q is the z interval occupied by the entire multilayered structure, A is the 
exchange coupling coefficient, and K is the anisotropy coefficient. The demagneti- 
zation tensor for a layer has only one element, D zz ; An is its value for an infinitely 
thin flat ellipsoid JTJj. In practice, one approximates Hi by the expression 

1 K 
Hi = H a + — [Ji t i +1 (m i+1 - mi) - Ji,i_i(m { - ra,_i)] - 2— ^e x x (m* x e x ) 



4irMi(mi ■ e z )e z , iel, 



(8) 
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where 

m = mi, m N+1 = m N . (9) 

The coupling coefficient J, which is related to A ( J = Ad~ 2 ), has the same value be- 
tween layers of the same material; similarly, the anisotropy coefficient K is constant 
within the same material, 

f Nh — 1, , , - .. v 

K h , i=l,... ,N h , 




Ji,i+i = < Jha, i = N h , K i = l »' N (10) 

The actual values of these material parameters depend on the temperature; K s <C 
Jf/j in all practical cases. 

2.2. Integration of the LLG equation. The LLG equation maintains a constant 
magnetization, so the only quantity that changes in the course of time is the direc- 
tion of the magnetic moment. We therefore begin by rewriting the LLG equation 
in terms of m. As the equation is entirely local to each layer, we drop the index i 
temporarily. We use the prime ' to denote differentiation with respect to time. 

Let H be the strength of the magnetic field, and let ft = H/H be the unit vector 
in the direction of H , 

H(t) = H(t)h(t). (11) 

Then the LLG equation is 

m' = —H [(m x ft,) + gm x (m x h)] . (12) 
We decompose the equation by means of the projection operators P and Q, 

Pu = (u ■ h)h, Qu = u - Pu = h x (u x h), u e R 3 . (13) 
Equation ([l2]) is equivalent to the two equations 

Pm' = -HP [(m x h) + gm x (m x h)} , (14) 
Qm' = —HQ [(m x h) + gm x (m x h)] . (15) 
Notice the identities 

P(m xh) = 0, P[m x (m x h)} = (m ■ Qm)h = -[1 - {Pm ■ Pm) 2 ]h, (16) 
Q(m X h) = —JQm. Q[m x (m x h)] = (m ■ h)Qm, (17) 
where J is the square root of the negative identity in R 2 , 

'-(i !)■'-(! -;)•*--'■ <*> 

Hence, we can recast Eqs. (|l4]) and ( p^| ) in the form 

Pm' = gH[l - (Pm ■ Pm) 2 ]h, (19) 
Qm' = i?[J - </(m ■ h)I]Qm. (20) 

Suppose that the direction of H does not change on an interval (t, t + At), 

h(s) = h(t), se{t,t + At). (21) 

Then Pm' = (Pm)' and Qm' = (Qm)' on (t,t + At), so Eqs. © and @ reduce 
to a coupled system of differential equations for the scalar u = (Pm ■ ft) in R and 
the vector v = Qm in R 2 , 

u' = gH(l-u 2 ) on (t,t + At), (22) 

v' = H(J-guI)v on (t,t + At). (23) 



HYSTERESIS IN LAYERED SPRING MAGNETS 



5 



From these equations we conclude that the critical states are u = 1, v = (m = h, 
magnetic moment parallel to the magnetic field) and u = — 1, v = (m = —ft-, 
magnetic moment antiparallcl to the magnetic field). The former is linearly stable, 
the latter unstable under infinitesimal perturbations. 

We now turn to the integration of Eqs. ( |2^ ) and (p3|). The former is independent 
of v and can be integrated immediately. If not only the direction, but also the 
magnitude of H is constant on (t, t + At), 

H(s)=H(t), se(t,t + At), (24) 

we find 

_ u(t) cosh(gH(t)( S - t)) + smh(gH(t)(s - t)) 
1 ' ~ cosh(gH(t)(s - t)) + U{t) smh(gH(t)(s - t)) ' fc 1 ' + j ' 1 j 



Next, we turn to Eq. (]2^). We replace the constant gH by — u 2 ) (from 

Eq. (p2|)) and use the identity —uu'/(l — u 2 ) = (ln(l — u 2 ) 1 / 2 )' to convert the 
equation into a differential equation for the vector w = (1 — u 2 ) _1 / 2 t>, 

w' = HJw on (t,t + At). (26) 
This equation can be integrated, 
u,( s ) = e 11 ^-^^^) 

= [cos(H(t)(s-t))I + sin(H(t)(s-t))J}w(t), se(t,t + At). (27) 
From the expression ( p5| ) we obtain 

\2\l/2 



(i _ , i(s) 2 )V 2 = (i-u(tyy 



(28) 



" (S) = cosh( g g ( t)( S -t)) + ^) S inh( g g(t)(^f)) " W ' se (^ + Ai )' ( M ) 



These results suggest the following choice of the integration scheme for Eq. ([L2| 
(m n ■ h n ) cosh(gH n At) + s'mh.(gH n At) 



cosh(gH n At) + (m n ■ h n ) s'mh(gH n At) 
cos(H n At)I + sin(H n At)J 



h„ 

h n x (m n x hn), (30) 



cosh(gH n At) + (m n ■ h n ) smh(gH n At) 
where m n+1 = m(t„+i), m„ = m(f„), /i„ = ^(i«), i?n = and At = 

The algorithm ( p0| ) is unconditionally stable for all values of At. Of course, 
the quality of the approximation suffers as At increases. However, the algorithm 
explicitly displays the relationship between the size of At and the local error in the 
time integration. The rate of precession of m around the polar axis is governed 
by H, the magnitude of the local effective field: in one time step, m precesses 
through an angle HAt. Therefore, by properly choosing At, we can resolve the 
fastest precessional motion in a given number of time steps per period. Since H 
varies over the course of a simulation, we have a natural and direct means to adjust 
the size of At to the current dynamical state, while maintaining the resolution of 
the precessional motion. 

Other algorithms for the numerical integration of the LLG equation have been 
proposed recently by Nigam Q and E and Wang Q . 
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2.3. Computing equilibrium configurations. The analysis in the preceding 
section suggests the following algorithm for finding equilibrium spin configurations. 
Starting from a given equilibrium state M = {Mi : i 6 /} at time to, one uses 
Eq. (||) to compute the magnetic field Hi in each layer at to. Having found Hi(to) 
for all i £ I, one advances in time to t± — to + At and uses Eqs. (§) and © to 
compute M at t\. If At is sufficiently small, M(t\) is a close approximation of the 
state of the system at time t\. One continues this process, finding approximations 
at successive times t n = t + rtAt, n. = 1,2,... , until equilibrium is reached. 



3. Numerical results. The algorithm of the preceding section has been used to 
study hysteresis phenomena in hard/soft bilayers that are driven by an applied 
field H a that is uniform, constant in time, and parallel to the planes of the atomic 
layers. The expression for the effective magnetic field, Eq. (|§|) , decomposes into an 
in-plane component, 

1 r 

[Ji,i+i{m i+ i - m l 



H l x e z = H a x e z 



Mi 
2 — m; 



J i ^i(m l - m;_i)] x e z 



el, 



(31) 



and an out-of-plane component, 



H, 



= [Ji,i+i(m i+ i 



mi) - Ji,i-i(mi - m,^)] • e z 



— 2 — m; 

Mi 



AnrMirrii ■ e. 



e I. 



(32) 



When the system is in an equilibrium state, the effective magnetic field is parallel 
(or antiparallel) to the magnetic spin; see Section 2.2. Hence, each Hi is a multiple 
of m^, and Eq. d32| ) reduces to a homogeneous system of linear algebraic equations 
for the set of scalars {m^ • e z : i G /}. In general, this system admits only the 
trivial solution, so the magnetic moments lie in the plane of the atomic layers. In 
the notation of Eq. (||), fa = for all i G I at equilibrium, and the only relevant 
variables are the in-plane angles {9i : i G /}. (Of course, the magnetic spin may 
have an out-of-plane component during the transient phase of the computation.) 

In the numerical simulations we focus on the in-plane angle of the magnetic spin 
at equilibrium and investigate its behavior as a function of the strength H a and the 
direction 8 a of the applied field, 



H a = H a h a , h a = (cos a , sin 



(33) 



The following computations refer to a bilayer configuration consisting of Nh = 115 
atomic layers of Sm-Co (a hard material) and N s = 100 atomic layers of Fe (a 
soft material). A different configuration is used in Section |3.3| , where we make a 
comparison with some magneto-optical measurements. Table [I] gives the values 
of the material parameters A, K, and M, as well as the values of the coupling 
coefficient J (J — Adr 2 , d = 2 A). In all cases, the damping coefficient g — 0.5. 

Table 1. Numerical values of the parameters. 





A (erg/cm) 


J (erg/cm 3 ) 


K (erg/cm 3 ) 


M (emu/cm 3 ) 


Fe 


2.8 ■ icr b 


7.0 ■ 10 M 


1.0 ■ 10 3 


1,700 


Interface 


1.8 • 1(T 6 


4.5 ■ 10 9 






Sm-Co 


1.2 • 10~ 6 


3.0 ■ 10 9 


5.0 ■ 10 7 


550 
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3.1. Rotational hysteresis. The case H a — 4800 oersteds is typical, at least for 
moderate values of H a (see Section 3.2). 

The simulations show that the equilibrium spin configurations for increasing 9 a 
(0 < 9 a < 2n) and decreasing 9 a (2n > 9 a > 0) are mirror images of each other. 
Figure [l] shows two sets of magnetic spin configurations at equilibrium for various 
values of 9 a , one set (left) as 9 a increases from to 2tt, the other set (right) as 9 a 
decreases from 2tt to 0. The heavy dots represent the endpoints of the magnetic spin 
(a unit vector) in each layer for various angles 9 a ; the values of 9 a , in degrees, are 
indicated near the top layer. (The dots merge into a solid line where the magnetic 
spins in adjacent layers are close.) 




Notice that the chirality ( "handedness" ) of the chain of magnetic spins changes 
from positive at 9 a = 301.5 to negative at 9 a = 301.6 degrees and from negative at 
9 a = 58.5 to positive at 9 a = 58.4 degrees. Figure || shows this change in a different 
way. 




-90 90 180 270 

In-plane angle of effective field, (degrees) 

FIGURE 2. In-plane angle 9i vs. i; H a — 4800 oersteds; (a) 9 a = 45, 
(b) 9 a = 90, (c) 9 a = 135, (d) 9 a = 180, (e) 9 a = 225, (f) 9 a = 270, 
( g ) 9 a = 301.5, (h) 9 a = 301.6, (i) 9 a = 315 degrees. 

Here, we have plotted the in-plane angle 9i against the layer index i for increasing 
values of 9 a . (The graphs for decreasing values of 9 a are obtained by symmetry.) 
First, the graph changes continuously (but not monotonically) as 9 a increases from 



8 



J. SAMUEL JIANG, HANS G. KAPER, AND GARY K. LEAF 



to 301.5 degrees, 9i increasing with i (positive chirality). Then it changes discon- 
tinously as 9 a increases to 301.6 degrees: 9i suddenly becomes decreasing instead 
of increasing with i. Finally it changes continuously again as 9 a increases further, 
9i decreasing with i (negative chirality), to return to the original graph (9i = for 
all i £ I) as 9 a reaches 360 degrees. In all cases, the spin is fixed along the easy 
axis (9i = 0) in most of the hard layers; it begins to deviate from the easy axis only 
as one approaches the interface (i — 115). The first derivative is discontinuous at 
the interface, and the tangent is vertical in the top layer (i — 215). 

The change in chirality is irreversible and induces rotational hysteresis. The in- 
plane angle of each spin vector traverses a different trajectory as the applied field 
rotates 360 degrees in the forward and backward direction. The hysteresis loop 
has the same shape, and particularly the same width, in all layers. Its vertical 
dimension contracts gradually as one descends through the soft layers, to disappear 
entirely in the hard layers somewhat below the interface; see Fig. |^. 



360 




90 180 270 360 

In-plane angle of applied field, 9 g (degrees) 

Figure 3. Rotational hysteresis: in-plane angle 6% vs. 9 a ; H a = 
4800 oersteds; i = 95, 115, 135, 155, 175, 195, 215. 



3.2. Two types of rotational hysteresis. When the strength of the applied field 
is varied, we observe different modes of irreversible behavior. We recall (Fig. ||) 
that, as 9 a increases from 0, the chirality changes discontinuously from positive 
to negative as the direction of the applied field deviates sufficiently from the easy 
axis. We denote the critical value of the angle 9 a by 9 C (9 C — 301.5 ... at H a = 
4800 oersteds). Figure f| shows the variation of 9 C with H a . 

As long as H a is sufficiently small, the magnetization process is reversible. At a 
first critical value of H ai marked H c \, the chirality of the chain of magnetic spins 
changes for the first time, and rotational hysteresis of the type discussed in the 
preceding section (with a basic period of 360 degrees) sets in. The width of the 
hysteresis loop, which is symmetric around 9 a = 7r, increases monotonically from 
at H a = H c i to some value less than 2-7T. 

At a second critical value of H a , marked H C 2, a sharp discontinuity occurs. 
The hysteresis loop narrows significantly and continues to narrow as H a increases 
beyond H C 2- The cause of this discontinuity becomes obvious in Fig. [jj, where we 
have plotted 9i against i; cf. Fig. |[ (The bottom 80 layers of hard material, where 
9i does not deviate noticeably from 0, are not included in this figure.) As H a 
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2 4 6 8 10 12 

Strength of applied field, H (kOersted) 

Figure 4. The critical angle C as a function of H a - 

reaches the value H C 2, the chain of spins has been stretched to its widest extent; 
it can no longer support the span in the top layer, stiffens suddenly, and becomes 
more like a rigid rod. The rod-like behavior is apparent from the increasing range 
where the chain is almost vertical. 




FIGURE 5. In-plane angle 0; vs. i; (a) H a = 1000, (b) H a = 2000, 
(c) H a = 3000, (d) H a = 5000, (e) H a = 6000, (f) H a = 7000, 
(g) H a = 8000, (h) H a = 9000, (i) H a = 10,000 oersteds. Right 
branches: 9 a just below C , left branches: 6 a just above C . 

The structural change in the chain of spins has some of the characteristics of a 
phase transition. For example, we observe a significant increase in the equilibration 
time (by two orders of magnitude) as a approaches C ; see Fig. |^. Also, the 
increasing size of the rigid domain near H C 2 is reminiscent of a diverging correlation 
length. 

At a third critical value of H a , marked H c a, another significant change occurs. 
The field has now become sufficiently strong to move the spins in both the soft 
and the hard materials. From here on, the chain of spin vectors changes over its 
entire length, maintaining its chirality. Figure ^ gives 0j vs. i for increasing values 
of 6 a . (The graphs for decreasing values of 9 a are obtained by symmetry.) The 
value H a = 10,400 oersteds is just above H c3 . This figure should be compared 
with Fig. for the standard case, H a = 4800 oersteds. The exact determination of 
H c3 is delicate; in our numerical simulations we found a slight rate dependence in 
the regime near H c3 . 
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Figure 6. Equilibration time near 6 C ; H a — 6797 oersteds, 9 C = 
323.0 degrees. 
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In-plane angle of effective field, 6 (degrees) 



Figure 7. In-plane angle 9 l vs. i; H a = 10,400 oersteds; (a) 9 a = 
16, (b) 8 a = 61, (c) 9 a = 106, (d) 9 a = 151, (e) 9 a - 195, (f) 6 a = 
196, (g) 6 a = 241, (h) 6 a = 286, (i) 6 a = 331, (j) 6 a = 375 degrees. 



Because the chain of spins behaves more like an elastic spring than a stiff rod, a 
new type of rotational hysteresis emerges, whose basic period can be any multiple of 
180 degrees. Figure || shows three graphs: one graph (c) is along the diagonal; the 
other two (a and b) are symmetric with respect to the diagonal. The outer graph 
(a) shows 9i for i = 85 (hard layer). The part below the diagonal is traversed in the 
upward direction as 9 a increases from 0; the part above the diagonal is traversed in 
the downward direction as 8 a decreases from 360 degrees. The spin is oriented in 
either the positive or the negative x direction. Transitions occur at 8 C and at every 
multiple of tt beyond 6 C . The center graph (c) shows &i for i = 215 (top layer). 
The orientation of this spin varies continuously with a and is perfectly reversible. 
Finally, the middle graph (b) shows 0t for % = 115 (at the interface). Here, the spin 
rotates continuously until it jumps. The jumps occur at 9 C and at every multiple 
of 7r beyond C . The graphs for the remaining layers fill the space between the ones 
drawn in the figure. The main point to observe is that the graphs for 9 a increasing 
always increase and stay below the diagonal, while those for 9 a decreasing always 
decrease and stay above the diagonal. Hence, chirality is preserved in both cases. 

When the direction of 9 a is reversed, 0j crosses the diagonal as soon as 9 a — 9 C is 
a multiple of 7r; after crossing, it remains on the part of the graph situated on the 
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90 180 270 360 

In-plane angle of applied field, 6 a (degrees) 

Figure 8. In-plane angle Q l vs. 9 a ; H a = 10,400 oersteds; (a) i = 
85, (b) i = 115, (c) % = 215. 

newly reached side of the diagonal. Because there is a gap between the graphs for 
8i in the interior layers and the diagonal, the orientation of the magnetic spin shows 
rotational hysteresis in all interior layers. This hysteresis is caused by a full-length 
transition of the chain of magnetic spins, rather than the partial-range transition 
that was responsible for the hysteresis below H c ^. 

3.3. Comparison with experiment. Quantities such as the magnetic moment 
are fundamental to describe the state of the system, but they are not directly 
measurable in an experiment. Measurable quantities are the torque density T and 
the magnetization angle a, 



T = H a d^2,M, 

iei 



j), a = tan 



-i J2iei MiS.YD.6i 
J2ieI M i C0S e i 



(34) 



Both T and a reflect the hysteretic behavior of the magnetic moments. Figure 
shows the torque density computed at H a — 4800 and 10,400 oersteds. 




90 180 

In-plane angle of applied field, 



270 

(degrees) 



Figure 9. Torque density; H a 
oersteds. 




90 180 270 360 

In-plane angle of applied field, 6 g (degrees) 

4800 (left) and 10,400 (right) 
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Experimental torque measurements at comparable values of H a show similarly 
shaped graphs, with extrema at approximately the same values of 9 a , but signifi- 
cantly narrower hysteresis loops Jl0[ . 

In Fig. we compare results for the magnetization angle with experimental 
data. The data were obtained by magneto-optical means for a bilayer consisting 
of Nh = 100 atomic layers of Sm-Co and N s = 250 atomic layers of Fe; the simu- 
lation curves also refer to this configuration p"T[ . The measurements were done at 
relatively low fields (H a = 360, 600, and 840 oersteds) and for a limited range of 
directions (6 a = : 10 : 230 degrees). 




90 180 270 360 

In-plane angle of applied field, 9 (degrees) 



FIGURE 10. Magnetization angle; (a) H a = 360 (o), (b) H a = 600 
(+), and (c) H a = 840 (x) oersteds. 



There is certainly qualitative agreement, but the simulations generally yield 
wider hysteresis loops than the experiments, and the discrepancy becomes greater 
as the field strength increases. This behavior can be explained by the fact that the 
model used in the simulations is a single-domain model, which does not allow for 
the important phenomenon of nucleation and motion of nanodomains. As a result, 
the demagnetization energy is seriously overestimated. In realistic simulations, one 
must use multidimensional models and allow for lateral inhomogeneities [ fTl") . 

4. Conclusions. In this article we have addressed an important issue in micro- 
magnetics: magnetization reversal in layered spring magnets. We have used a one- 
dimensional model of a film consisting of atomic layers of a soft material on top 
of atomic layers of a hard material with strong coupling at the interface, assuming 
no variation in the lateral directions. The state of such a system is described by a 
chain of magnetic spin vectors. Each spin vector behaves like a spinning top driven 
by the local magnetic field and subject to damping. The dynamics are described 
by a system of LLG equations, Eq. (g), coupled with a variational equation for the 
magnetic field, Eq. (||). 

We have presented an integration procedure that maintains the invariance of 
the magnetization (the magnitude of the magnetization vector) and proposed an 
algorithm for finding the equilibrium state of the system. 
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We have applied the algorithm to simulate magnetization reversal in layered 
spring magnets. The results show that a layered spring magnet exhibits rotational 
hysteresis with a basic period of 360 degrees at moderately strong fields and ro- 
tational hysteresis with a basic period of 180 degrees at strong fields. The former 
type of hysteresis is induced by a partial-length transition of the chain of magnetic 
spins; the transition occurs only in the soft material and causes a change of chiral- 
ity. The hysteresis in strong fields is induced by a full-length transition of the chain 
of spins in both the hard and the soft layers; it is much weaker than the rotational 
hysteresis at moderately strong fields and can cover any period that is a multiple 
of the basic period. 

The numerical results for the torque and magnetization angle agree qualitatively 
with the experimental data but differ at the quantitative level. In particular, the 
one-dimensional model seriously overestimates the demagnetization energy since it 
does not allow for the nuclcation and motion of nanodomains. In realistic simula- 
tions, lateral inhomogeneities must be taken into account. 
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